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Abstract 

Despite the seminal connection between classical multiply-periodic motion 
and Heisenberg matrix mechanics and the massive amount of work done on the 
associated problem of semiclassical (EBK) quantization of bound states, we 
show, that there are, nevertheless, a number of previously unexploited aspects 
of this relationship that bear on the quantum-classical correspondence. In 
particular, we emphasize a quantum variational principle that implies the 
classical variational principle for invariant tori. We also expose the more 
indirect connection between commutation relations and quantization of action 
variables. In the special case of a one-dimensional system a new and succinct 
algebraic derivation of the WKB quantization rule for bound states is given. 
With the help of several standard models with one or two degrees of freedom, 
we then illustrate how the methods of Heisenberg matrix mechanics described 
in this paper may be used to obtain quantum solutions with a modest increase 
in effort compared to semiclassical calculations. We also describe and apply 
a method for obtaining leading quantum corrections to EBK results. Finally, 
we suggest several new or modified applications of EBK quantization. 

21.60.-n,21.60.Ev,21.60.Jz,03.65.Ca 
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I. INTRODUCTION 



Though applications of great interest (and increasing complexity) continue to be de- 
veloped, e. g., fl|||, the theory of the semi-classical quantization of invariant tori by the 
application of EBK quantization conditions appears to be a closed (or at least qui- 

escent) sector in the study of the relationship between the quantum mechanics and the 
classical mechanics of non-separable systems. (For exceptions, see the work of Jaffe and 
collaborators and the even more recent work of Morehead f9j.) Focusing on systems 

with two degrees of freedom, we can divide the many methods that have been developed 
and applied to this subject into two main sub-categories, those based on the solution of 
Hamilton's equations as an initial- value problem, i. e., the calculation of trajectories, and 
those based on algebraic methods involving trajectories indirectly or not at all. Within this 
subdivision, referred to as Ml and M2, respectively, we may distinguish principally: 

Mia) Independent actions are computed from closed curves, signaling invariant tori, gener- 
ated by the intersection of trajectories on two independent surfaces of section |T(]|-|T2"|] . 

M2a) If invariant tori exist, the dynamical variables can be represented as multiply-periodic 
functions of the angle variables, with Fourier coefficients that depend only on the 
frequencies, or equivalently, actions. Equations of motion for the Fourier coefficients 
are obtained in either Hamiltonian or Lagrangian form, and for each there is an asso- 
ciated variational principle. The equations of motion are solved either perturbatively 
or non-perturbatively and the quantized actions and energy calculated in terms of 
the known Fourier coefficients |T3|-pT|. This is the approach of paramount interest 
to us. 

Mlb) Consider the system Hamiltonian to be a sum of an integrable part and a pertur- 
bation. The latter is turned on adiabatically over a time T. From the assumption 
that the actions computed initially for the unperturbed system are approximate adi- 



abatic invariants, it is possible to obtain quantum energies [p2|,p3| (sometimes even 



for values of the coupling strength at which associated invariant tori no longer exist). 
M2b) The Hamilton- Jacobi equation is solved iteratively, using a Fourier series expansion, 



thus providing the generating function for the appropriate action-angle variables |24 
(In contrast the Fourier series constructed in M2a are to be understood as the explicit 
equations of transformation from the original dynamical variables to the action-angle 
set.) 



Mlc) It can be shown that a finite set of adjacent non-quantized trajectories with 

the same total energy can be used to calculate accurate values of the actions for 
one of them. Quantized energies and associated actions are computed by linear 
extrapolation. In this method it is necessary to propagate trajectories until they 
almost close on themselves. 



M2c) Following the ideas of Birkhoff and Gustavson p7j , several groups pq-pT| have car 



ried out increasingly ambitious programs for transforming a given Hamiltonian to 
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normal form by a succession of canonical transformations. The resulting Hamiltonian 
is quantized and energy values obtained. 



Mid) Trajectories propagated over a sufficiently large multiple of the elementary periods 
of a multiply-periodic orbit can be Fourier transformed to yield the Fourier repre- 



sentation of the fundamental dynamical variables. Applying Percival's formulas |I7 
the actions are computed. Initial conditions are varied until quantized orbits are 

found mm. 



For further discussion and a more exhaustive list of references, we refer the reader to 

Refs. mus 

Despite all this effort, there remains a gap in the study of the relationship between 
quantum mechanics and the theory of invariant tori, especially as it relates to the work that 
utilizes the description of invariant tori in terms of Fourier series. This assertion is grounded 
on the observation that each form of quantum mechanics is associated in a natural way with 
a corresponding version of classical mechanics. Thus in the same sense that Schrodinger 
wave mechanics is naturally associated with the Hamilton- Jacobi equation, and the classical 
limit of the path-integral method is Hamilton's variational principle, it is equally true that 
the classical limit of Heisenberg matrix mechanics for bound systems are the equations of 
motion for multiply-periodic systems. Though some of the work cited is based, in the sense of 
the correspondence principle, on this passage between quantum and classical descriptions, 
and the passage itself is in some respects well known (see for instance), the main 

thesis of this paper is that this route for passing from the quantum to classical theory has 
not been fully explored. A possible explanation is that although the basic correspondence 
is well-known, it is hardly known outside the circle of the authors and former associates 
that matrix mechanics can be derived from a novel variational principle [j36|-39l called the 



trace variational principle (see also [|40|| ) and that the classical limit of this principle is the 
variational principle for invariant tori 

Guided by this relationship, the main purpose of this work is to study the transition 
from the quantum to the classical domain more thoroughly than heretofore for both the 
dynamics and the kinematics (quantization conditions). This has led to some results that 
we find it difficult to believe are not known, but for which we do not have an independent 
reference. In Sec. II, we supply a brief but self-contained account of our version of matrix 
mechanics, with emphasis on the variational basis. The passage to the semi-classical limit 
is then carried out in Sees. Ill and IV, where we show that the limit of the quantum 
variational principle is the variational principle for invariant tori. We find, however, that 
the relationship between the commutation relations and the EBK quantization conditions is 
more indirect, the former corresponding in the classical limit to the Poisson bracket relations 



and the latter to the Lagrange bracket relations In Sec. V, we study the commutation 
relations per se, the most tangible result being a new derivation of the WKB quantization 
condition, applicable to one- dimensional systems. Section VI is devoted to some illustrative 
numerical studies. We describe in turn and then apply algorithms for carrying out the 
semiclassical calculations, the associated matrix quantum calculations, and a method for 
calculating directly quantum corrections to the semiclassical result. The need for the latter 
as a separate approach arises from the fact that the matrix quantum calculations are designed 
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to give better results than the semiclassical one for low-lying states but do not go over in any 
limit to the EBK calculation. Finally, Sec. VII contains several suggestions for new ways to 
use the semiclassical quantization scheme. In a final section, we make some proposals for 
further work. Preliminary accounts of the main new theoretical results of this paper can be 
found in HQ. 



II. MATRIX MECHANICS 

A. Variational principles and equations of motion 

We begin with a brief description of a variational principle, the so-called trace variational 
principle, from which one can derive Heisenberg's form of quantum mechanics. Though a 
version of this principle was suggested more than three decades ago f3~6|l , and subsequently 



several publications have been devoted to its exposition and further development |37H39j, 
it appears to be largely unknown by the community at large. Except for one brief allusion, 



2jJ , it has not been applied to the problem at hand. 

Most of the important elements are already present for a system with one degree of 
freedom, and we therefore focus attention on the Hamiltonian 

H= l -p 2 + V(x), (2.1) 

with equations of motion 

[x,H]=ip, (2.2) 



[ip, H] = dV/dx = V, (2.3) 
derived by utilization of the commutation relation (U — 1) 

[x,p]=i. (2.4) 



In practice we are usually concerned with the matrix elements of ( |2.2j )-( |2~4l ) in the rep- 
resentation in which H is diagonal, with eigenvalues E n . Namely, 



J mni 



\E n E m )x mn ipn 
{E n E m )ip mn = (V ) m?1 , (2-5) 

and 

[x,p)nm=i$nm, (2-6) 

where x nm = (n\x\m). We shall feel free to use both notations interchangeably. 



In early work [ J3| we have shown, for polynomial potentials, how the energy differences 
and the matrix elements x mn , p mn can be obtained from Eqs. ( [2.5|) and (|2.6|). The primacy 
of these elements can be seen by using the completeness relation for the evaluation of matrix 
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elements of a product. The eigenvalues themselves can be found by the direct evaluation of 
the expectation values 

E n = H nn = (n\H\n) = ]T ^|p nn f + (n\V(x)\n), (2.7) 

n' 

where the application of completeness is illustrated in the kinetic energy term. Ultimately, 
we shall be concerned with extending the previous algorithms to the multidimensional case. 

Before proceeding to the discussion of a variational formulation, we add a few remarks 
about the implementation of the above formalism. Equations ( |2.5|) and (|2.6| ) are, to start 
with, an infinite set of sum rules that must be satisfied by the exact eigenstates of the 
Hamiltonian. Starting at any point in the spectrum, sums spread without cutoff as far as the 
configuration space will allow. To obtain closure, we must make two kinds of approximation. 
The first is that the matrix elements are rapidly decreasing functions of k = |n — m\, so 
that all matrix elements with k > k max can be set to zero. The second is that the matrix 
elements are, for sufficiently large n and m, slowly varying functions of n + m, so that for 
some values sufficiently far from the center of interest in a given calculation we can set 
(n + r\x\m + r) = (n\x\m) for |r| <C \( m + n )- With these assumptions Eqs. fl2.5| ) and 
( p.6| ) reduce to a finite set of equations, where in addition to the retained matrix elements 
of x and p, the additional unknowns are a set of elementary energy differences from which 
all other energy differences can be composed. To determine this collection of variables, it 
suffices to utilize all the available equations of motion, but only the diagonal elements of 
the commutator. This is consistent with the result, proved in the next subsection, that the 
off-diagonal elements of the commutator are a consequence of the equations of motion. 

A natural question to ask is whether Eqs. (|2.5| ) can be derived from a variational prin- 
ciple? Here we wish to treat the matrix elements of x and p as variables in the variational 
statement 5E n = 5H nn = 0. There are, however, two obstacles to such an endeavor: (i) The 
matrix elements are not all independent, (ii) The same matrix elements appear in different 
energy functionals. Thus p nn i occurs both in H nn and in H n i n i. For which is it to be a 
variational parameter? A solution to the second problem posed is to form an average of the 
stationary functionals. It turns out that in order to derive the equations of motion as given 
above, it is necessary to choose the most symmetrical possible average, namely, the trace. 
Thus we require 

8j2H n n = STrH = 0. (2.8) 

n 

A solution to the first problem is to impose all the possible kinematical constraints, namely 

5[x,p) nn , = 0. (2.9) 

Multiplying (|2.9|) by a Lagrange multiplier matrix (— i)A n i n (A is Hermitian), we add the 
result to (|2.8| ) and are thus led to a master variational principle 

= 5F 
= 5Tr{H -iA[x,p]} 
= 5Tr{H -ip[A,x}} 

= 5Ti{H + ix[A,p]} . (2.10) 
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The several forms are equivalent because of the assumed cyclic invariance of the trace. (This 
is certainly unobjectionable in practice where the trace is taken over a finite dimensional 
vector space.) 

Carrying out the unconstrained variation of (|2.10|) with respect to the matrix elements 
x n > n and Pn'n, keeping A fixed, and using the explicit form ( |2.1| ) of H we obtain the equations 

Pn'n = A]nn', ( 2 -H) 

(V') nn , = i\p,A} nn ,. (2.12) 

Because of the invariance of the trace with respect to choice of basis, the representation |n) 
is, at this point, arbitrary. The most convenient immediate choice is the one in which the 
Hermitian operator A is diagonal. By comparing with the known equations of motion, we 
then identify A as the Hamiltonian. 

The derivation of the equations of motion does not exhaust the consequences that can 
be drawn from the trace variational principle. We shall now demonstrate from this principle 
that the vanishing of the off-diagonal matrix elements of the canonical commutator is a 
consequence of the equations of motion, leaving only the diagonal elements as independent 
kinematical conditions. This result is consistent with the "empirical" finding above that 
a dynamical scheme is, in fact, completely determined by adjoining this one-dimensional 
(diagonal) array of kinematical constraints to the equations of motion. 

To derive the off-diagonal elements of Q2.6|) , we make use of the invariance of the trace 
under an infinitesimal change of basis. In the new basis, the Hamiltonian will not be diagonal, 
in general, and thus we must allow for a change in the Lagrange multiplier matrix. We 
calculate 

= 5F = Tr{-i(5H)[x,p]}, (2.13) 

since all other contributions vanish in consequence of the equations of motion. If we express 
the infinitesimal change of basis in the standard form 

S\ n >= -ieQ\n >, (2.14) 

where e is infinitesimal and G is Hermitian, we recognize that in a variation about the energy 
diagonal representation 

6 < n\H\n >= ie < n\[Q, H]\n >= 0. (2.15) 

On the other hand non-diagonal elements 

5 < n\H\m >=< n\5H\m >= 5H nm (2.16) 

can be assigned arbitrary infinitesimal values consistent with Hermiticity. From this and 
( [2.131 ) we conclude that the off-diagonal elements of [x, p] vanish. 

One additional result of great importance is that a solution of the dynamical scheme 
proposed above guarantees that the Hamiltonian is diagonal, H n , n ' = 0, n' . This result 
is derived in the Appendix. 
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In place of Hamilton's equations ( |2.5| ) it is often more convenient to consider Newton's 
equation 

This equation may be derived from its own variational principle by substituting the first 
of Eqs. ( p.5|) into the previous functional F. The result is a new functional G that can be 
written in the alternative forms 

G = Tt{H - H[x, [H, x}]} 

= Tr{~[x,H][H,x) + V(x)}. (2.18) 

In the second form, we recognize that G is the negative of the trace of the Lagrangian. To 
obtain Newton's (or Lagrange's) equations one varies G with respect to the matrix elements 
of x, keeping the matrix elements of the Hamiltonian fixed. From the structure of fl2.18| ), 
this means keeping the energy differences fixed. This formulation of the quantum theory 
is completed by adjoining a form of the canonical commutator from which the momentum 
operator has also been eliminated, namely 

= J2(2Ei-E m -E n )x ml x ln . (2.19) 
l 

It should also be mentioned that the vanishing of the off-diagonal matrix element of the 
commutator can equally well be proved by paraphrasing for the functional G the argument 
presented for the functional F. 

Since our main interest in this paper is in autonomous systems with at least two degrees 
of freedom, we must now describe how the previous considerations are modified by this 
generalization. We therefore consider a system with N coordinates x = (x±, ...,xn) and a 
Hamiltonian of the form 

H = Y.\pi+nx). (2-20) 

i 

The functional F from which Hamilton's equations are derived takes the form 

F = Tr{H-iA y £[x i ,p i ]}, (2.21) 

i 

where A is once again identified as the Hamiltonian. Lagrange's equations are derived from 
the functional G, where 

G = Tr{-~£[^, + V(x)}. (2.22) 

Before providing any further details, we must mention the problem of labeling of the 
eigenvalues and eigenstates of H. In our earlier work f21| , we "naturally" assumed that 
the labeling could be done by a choice of N integers n = (ni, njy) of which we could 
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keep track as we tuned one or more coupling parameters, starting from values for which the 
problem was integrable. The same assumption was discussed rather more thoroughly by 
Percival [p^-p!6| who emphasized that the validity of this assumption is coterminous, in the 
semiclassical limit, with the existence of invariant tori. 

With this understanding, it is not necessary to write the equations of motion a second 
time, but only to remember that there is now one equation for each value of i and in each of 
these equations to replace the integer n by the corresponding vector n. There does remain 
one question to be addressed that will be of some importance to us later. This is the 
question of whether we can deduce from the variational principle the separate vanishing, for 
each value of i, of the off-diagonal elements of the commutators This follows from 

the fact that the rows and columns of the matrix A are each labeled by an N- dimensional 
vector, provided the set of non-vanishing matrix elements of the product XiPi is disjoint 
from the corresponding set for any other choice of coordinate index. This will be true for 
any model that we study. We have been somewhat cavalier in the present discussion, but 
the questions that we have slighted, in particular, why other elements of the algebra do 
not appear in the constrained variational principle, will be considered in more detail in the 
discussion that now follows. 



B. Commutation relations and equations of motion 

The canonical commutation relations must be constants of the motion; this means that 



their time derivatives should vanish ||45|| . We show this for the class of Hamiltonians under 
study. Consider first 

= 0. (2.23) 

This calculation shows that in the energy- diagonal representation (assuming no degeneracy) 
the commutator of two different components of momenta has no off-diagonal matrix ele- 
ments. If we choose the momenta to be imaginary Hermitian operators, it follows that the 
diagonal elements of the commutator also vanish. Thus it is clear that the commutators 
[Pi>Pj] = rnay be omitted from the dynamical scheme. 

Utilizing the previous result that the off-diagonal elements of the commutator 
vanish, we calculate 

Since V is to a large extent arbitrary, we may safely conclude that the coordinates all 
commute with one another. 

Combining the previous two results, we next verify that 

d 

^fc'Pil = \P^Po] + -fa] = °- ( 2 - 25 ) 
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Finally we check compatibility by the calculation 

— [x h Xj] = [xi,pj] + \pi,xj] = 0. (2.26) 

We have thus shown that the commutation relations are compatible with the equations 
of motion for the class of Hamiltonians under consideration. For the practical problem 
of constructing a calculus based on Heisenberg matrix mechanics, the consequence of our 
deliberations is that at most only the elements diagonal in the energy representation of the 
commutators of a coordinate and the corresponding momentum can enter, if we use all the 
available equations of motion. 



III. SEMICLASSICAL LIMIT: MATHEMATICAL PRELIMINARIES 

The purpose of this section and the one to follow is to show that the semiclassical 
theory of invariant tori is the "natural" limit of the quantum theory of the previous section. 
We shall first collect in the form of lemmas some of the mathematical statements that we 
need. In the following we shall use the notation Lim (n|0|n) to signify the leading term 
in the semiclassical limit of the designated matrix element. Here O is generally a product 
of elementary operators. We consider first the one- dimensional case and define for a real 
Hermitian operator, A, 

A k {n) = (n- h\A\n + h) = A_ fc (n). (3.1) 

Notice that this quantity is an analytic continuation of a nearby physical matrix element 
and is our definition of the semiclassical limit of the matrix element (n\A\n + k) . As shown 
in a previous work [20], this choice can be used to provide a completely algebraic basis for 
the standard WKB quantization rule for a vibrational degree of freedom. 

We note that with the definition ( |3.1|) , we have for the first two terms of a Taylor 
expansion, 

, / s 1 , dA k (n) 
(n\A\n ±k)* A k (n) ± - k—^ 1 . (3.2) 

Next, with the array of amplitudes A k (n), for fixed n and varying k we associate a formal 
Fourier series 

oo 

A(n,9)= A k (n)exp(ik6), (3.3) 

k=— oo 

though in practice we shall always deal with severely restricted sums. Below, we shall then 
make extensive use of the average 



(ABC.) = (2tt)- 1 / deA{n,6)B{n,9)C{n,e)... , (3.4) 
which is just the constant term in the Fourier series of the product. 
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With the above preliminaries, we are prepared to state and prove a series of elementary 
propositions. 
Lemma 1: 

Lim (n\AB\n) = (AB). (3.5) 

For the proof we write 

Lim (n\AB\n) = Lim ^ [(n|A|n + k)(n + k\B\n) 

k>0 

+{n\A\n-k){n-k\B\n)]. (3.6) 

Here we have assumed that n is sufficiently large that the sum over k can be extended far 
enough in both directions, to k max , to obtain numerical convergence. This establishes the 
limited and non-rigorous sense in which the word proof is to be understood both here and 
below. 

If we now apply (|3.2|) in sequence first to the matrix elements of A, for example, and 
subsequently to the matrix elements of B, we find 

Lim (n\AB\n) = ^ A k (n)B_ k (n) [1 + 0(n~ 2 )] 

all A; 

= (AB)[l+0(n-% (3.7) 

The error estimate arises from the assumption that a derivative with respect to n is of 
relative order (1/n). That the error is of second order is a consequence of our choice of 
definition ([3j]). The sum on the right hand side of ([3. 71) has the value required by the 
lemma. 

Lemma la: With the same assumptions as before, Lemma 1 can be extended to a product 
of more than two factors, 

Lim (n\ABC...\n) = (ABC.) = ]T 5 kl+k2+ki+ ... , A kl (n)B k2 (n)C k3 (n)... . (3.8) 

The proof is carried out by ordering the various upward-going, downward-going and mixed 
contributions to the multiple sum so that Eq. (|3.2|) can be applied. 

Lemma lb: The previous lemmas can be extended to the mult i- dimensional case. Ex- 
tending the boldface notation now to designate, in addition to the quantum numbers n, also 
the integer vector k for the components of a multiple Fourier series and the vector for an 
array of angle variables, we introduce a formal multiple Fourier series, 

A(n,6) = ^A k (n)exp(zk- 6). (3.9) 

k 

The lemma then applies to the average 

(AB) = (2ny N J d0A(n,0)B(n,0) 

= £>k(n)B_ k (n) (3.10) 

allk 
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and to corresponding multiple products. 
Lemma 2: 



Lim(n\0\n + k) = Oy.{n) 

= W 1 J d6exp(-ik6)0(n,6). (3.11) 

Here is a product of two or more elementary operators, since for a single operator, the 
previous statement is only a combination of the definitions ( ^.1[ ) and ( |3.3| ). The same 
equation applies in boldface notation. We shall not actually need this lemma, since its 
application would be to obtaining the semiclassical limit of the equations of motion directly. 
Our procedure, however, will be to obtain that limit for the variational principle and then 
derive the equations of motion from the latter. 

We do require expressions for limits of commutators. We consider the usual Poisson 
bracket, 

[AB]pB = ^ [ ^ — &l ^ — mT ] - (3 - 12) 



We can now state 
Lemma 3: 



Lim{n\[A,B]\n + k) = i{[A, B] PB ) k , (3.13) 



where the notion on the right hand is understood as the hth Fourier component of the 
Poisson bracket ( |3.12|) . This lemma can be extended to multiple commutators. In practice, 



we shall only require the double commutator, for which we find 
Lemma 3a: 

Lim(n|[A,[S,C]]|n + k) = (i) 2 ([A, [B, C] PB } PB ) k . (3.14) 

We shall utilize these results to develop the dynamics and kinematics of a semiclassical 
quantization scheme. Before turning to this task, we note an important consequence of 
the results in this section. Whereas in quantities such as (n\ABC\n), different base values 
of n will occur, it is a consequence of the elementary reasoning applied above that in the 
expression Lim(n\ABC\n) only the reference value of n occurs. Thus in the semiclassical 
limit a trace decomposes into a sum of independent terms. 



IV. SEMICLASSICAL LIMIT: DYNAMICAL AND KINEMATIC AL SCHEME 

A. EBK Scheme 

We now apply the reasoning and results of the previous section to obtain the semiclassical 
limit of the (Lagrangian) form of the trace variational principle. We are interested in a 
value of the vector n for which the corresponding energy and eigenstate may be described 
accurately by the semiclassical approximation. We should then form the trace over a space of 
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states extending in both directions in choice of n compared to the reference state. From the 
results of the previous section, however, it follows that for the purpose of deriving equations 
of motion we may suppress the trace, because in the limit considered, as we have already 
pointed out and in contrast to the quantum case, the states decouple. It thus suffices to 
focus attention on a given state. (At the end of this section, however, we shall have occasion 
to restore the trace.) 

We consider first the kinetic energy T. If E(n) is the energy of the state identified by n, 
then here we shall encounter the correspondence principle in the form 

Ui(n) = E(m, ...rii + -, ...) - E(ni, ...m - -, ...) 

(dE(n)/dm), (4.1) 

i. e., the classical frequencies u>(n) = ...,a>jv) approach the quantum energy differences. 
We then find by a perfectly straightforward examination that 

Lim2(n|T|n) = LimJ](n|[x i , H][H, Xi]\n) 

i 

= LimJ2[E(n + k) - E(n)) 2 {n\xi\n + k)(n + k\xi\n) 

i,k 

= £(k ■ w) 2 x i , k (n)x i ,_ k (n) = (2T). (4.2) 

We note that the same result may be found by applying Lemma 3 for the classical limit of a 
Fourier component of a commutator, provided we recognize that the classical Hamiltonian 
is a function only of the action variables n. This follows from the fact that the formal 
Fourier series introduced to represent the solution of the classical problem defines a canonical 
transformation to the correct action angle variables (see below). 

For the quantity just evaluated, as well as for the potential energy, 

Lim (n|V|n) = (V), (4.3) 



the relative error is of order (1/n 2 ), provided we use the definition (|3.1|) . 
If we combine Eqs. ( |4.2p and ( |4.3[ ) we have 

Lim (n|L|n) = (L). (4.4) 

Since this expression is the semiclassical limit of a quantity stationary with respect to varia- 
tions of matrix elements of the coordinates, keeping energy differences fixed, we expect ( |1.4| ) 
to be stationary with respect to variations of the associated Fourier components, keeping 
the frequency components fixed. Conversely, requiring the stationary property, 

6{L) = 0, (4.5) 



Sx 



i,-k 



keeping the frequencies fixed, yields Lagrange's equations in Fourier component form, 
namely, 



12 



(k.u,) 2 z l)k (n) = |^. (4.6) 



Considered from the purely classical point of view, the solution of Eqs. (|4.6|) may be 
undertaken from several standpoints. The choice of real Fourier coefficients, adhered to 
throughout the present development, has already determined half of the initial conditions, 
namely, all components of the velocity vanish initially. If we then specify iV fundamental 
Fourier components of the coordinates, we thereby define a scheme, provided a solution 
exists, that determines the remaining Fourier components as well as the N frequencies. An 
alternative scheme, more integral to the variational principle, is to specify iV (incommen- 
surable) frequencies and calculate the Fourier amplitudes from the equations of motion. In 
this paper, however, our interest will be in adjoining iV quantization conditions to the equa- 
tions of motion. The resulting set, if it has a solution, must then determine the Fourier 
amplitudes and the frequencies. 

We apply the EBK quantization condition (restoring H for the instant) 

h = (th + -Oi)h, (4.7) 

where aj, the Maslov index, which has the value two for an uncoupled vibrational degree 
of freedom and for the models studied in Sec. VII. The action is given by the equivalent 
expressions 
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= ^2ki(k- u;)x k x_ k 

k 

= <p-J-*>, (4.8) 

where, in the last form, p is the momentum vector of the system. This completes the formal 
statement of the semiclassical quantization for invariant tori. 



B. Connection between EBK quantization and commutation relations 

We have, in fact, given two quantization schemes in this paper, the completely quan- 
tum one consisting of Heisenberg's equations of motion adjoined to diagonal elements of the 
canonical commutation relations and the semiclassical one consisting of the classical equa- 
tions of motion in Fourier component form adjoined to EBK quantization conditions. But 
a part of the logical story is clearly missing. On the one hand, we have derived the classical 
equations of motion from the quantum ones using correspondence principle arguments mod- 
ified only minimally from historical forms. On the other hand, although EBK quantization 
can also be derived, up to the value of the Maslov index, from the correspondence principle, 
what has not been made clear is the relationship between the canonical commutation rela- 
tions and the EBK quantization conditions. Though it is hard to believe that this question 
has not been addressed at some point in the history of quantum mechanics, we have never 
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encountered such a discussion in the literature. Therefore, for the sake of completeness, we 
include it at this juncture. 

Consider first the diagonal matrix element of the commutation relation 

Lim (n\[xi, [H,Xi]]\n) = -zLim (n|[xj,pj]|n) 

__ ,^rdandpi_ _ dxj dpi 

y Do, <)ii, dnjde^' 

d 

= V] tj kj (k • u;)xj ik £i._ k . (4.9) 

The basic result exhibited here, that follows from Lemma 3, is that the limit of the diagonal 
element of the commutator is the time average of the corresponding Poisson bracket (PB). 
On the other hand, we observe that the time average of the fundamental Lagrange bracket 
{9i,rii} has the value 

» i,k 

= = L < 410) 

Thus the EBK quantization condition is related to the Lagrange bracket (LB) rather than 
to the PB. 

It is, of course, well-known [51] that one set of fundamental brackets implies the other. 
But for this purpose we must consider the full brackets rather than just their time averages 
over the torus. We illustrate the argument for one degree of freedom. We thus wish to show 
that the formal Fourier series 

x(n,8) = Xgjn) exp(iq6), 
i 

p(n, 8) — }Jqux g (n) exp(iq9), (4-11) 
i 

satisfy the PB condition 

[x,p}pb = [x,p]e,n = 1- (4.12) 

provided, of course, that Eqs. (]4.11|) satisfy the equations of motion. In that case the Fourier 
series represent the canonical transformation between the given canonical coordinates x,p 
and the exact angle-action variables 6, n. By working out the Fourier series for the PB, 
requiring that the constant term be unity and that the other Fourier components vanish, we 
find, with the help of a linear transformation of the summation indices, that 

Id Id 
5 q ,o = EP 2 - l^ x k+^g^( UJX ~k+i q ) + {k~ ^qfujx_ k+ i q —x k+ i q . (4.13) 
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In the sum, k takes on both integral and half integral values. For q = this expression 
reduces to the one- dimensional form of (|4.9f) . The vanishing of the remaining Fourier com- 
ponents must be a consequence, we suspect, of the vanishing of the off diagonal elements of 
the commutator, i. e. , it must be the classical limit of this property. This can be verified 
directly, using Lemma 3a. As we emphasized at the beginning of this discussion, the equa- 
tions of motion must be involved, as indeed they are, in the vanishing of these off-diagonal 
elements. 

Using the results derived in Sec. IIB, we can extend the detailed considerations just 
given to the mult i- dimensional problem. This is in fact the interesting case, since in the 
one-dimensional case the PB and LB are indistinguishable, a fact that will be used later to 
derive the WKB quantization condition from the commutation relation. The correct value 
of the Maslov indices must somehow also be implied in the multi-dimensional derivation of 
the EBK quantization conditions. The only derivation we can supply at the moment is one 
based on starting with a separable system and assuming adiabatic invariance as one turns 
up the coupling. 

It may be of some interest to derive the q 7^ part of ( 4.13| ) directly from the classical 



variational principle. For this purpose, we find it necessary to retain the full structure of 
the quantum variational principle, so that in the classical limit we have not only an average 
over a given torus, but in addition, from the trace operation, an integral over all tori. In 
other words, we consider an average over phase space and indicate it by a double bracket 
notation, e. g., 

((L))= [ dnd6L(6,n), (4.14) 



where the phase average of L has the form, 

< L > = ? 2 «§> 2 <f > ! > " <">■ < 4 - 15 > 

We now subject the system to an arbitrary infinitesimal canonical transformation about 
the exact solution. In consequence of the equations of motion, the only non-vanishing 
contribution comes from the first term of (|4.15 ) and has the form 



sm) = -« ( §)*f*f» 

Here the second line has been obtained by an integration by parts with respect to n with 
boundary terms dropped. This requirement is clearly the classical analogue of the invari- 
ance of the trace, as is the requirement that (|4.16|) vanish. Finally if we assume that SH 
is expanded in a Fourier series with no constant term and otherwise arbitrary Fourier co- 
efficients, we are forced to the conclusion that all Fourier components of the expression in 
square brackets, other than the constant term, vanish. It is in fact easy to see that this 
conclusion coincides with the corresponding statement contained in ( 4.13Q . 



It is, finally, important to refer to the classical limit of the property that the quantum 
theory is formulated in the representation in which H is diagonal. The limit of this property 
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is precisely the condition that when the canonical transformation to the exact action-angle 
variables has been carried out, it is signaled by the vanishing of the Fourier series for the 
classical H, other than the constant term. In applications, we shall utilize this requirement 
in both its quantum and classical aspects as a test of the convergence of our solutions. For 
example, in the classical limit, the difference 



should vanish. 

We summarize in Table I the analogy that has now been fully established between Heisen- 
berg matrix mechanics in the energy- diagonal representation and the theory of invariant tori. 
For the sake of simplicity, the notation appropriate to one degree of freedom has been uti- 
lized. All entries have been described in the preceding text. Note that we have included the 
Hamiltonian rather than the Lagrangian form of the variational principle. 




(4.17) 
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TABLES 



TABLE I. The concepts of the theory of invariant tori as the classical limit of the concepts 
associated with Heisenberg matrix mechanics in the energy-diagonal representation. 



Quantum 



Semi- Classical 



1. 
2. 
3. 

4. 
5. 



9. 



Matrix elements of x,p 
(n\x\m), (n\p\m) 
Trace variational principle 
STr(H -ip[H,x\) = 
Matrix element of EOM 

{E n+ k 



En) %n+k,n 



(dV\ 

V dx Jn+k,k 

Commutation relations 
(n|[x,p]|n) = i 
Energy is diagonal element of H 
E n = (n\H\n) 
Quantum frequency 

<^n+fc,n = E n+k ~~ E n 

Operator 
A{x,p) 
Matrix elements 
(n|^4|n + k) 
Energy is diagonal 
J2kM n \ H \ n + k )( n + k\H\n) = 



Fourier components of x(I,9),p(I,0) 

Xk{I), Pk{I) 
Constrained variational principle 

6({H- U p%))=0 
Fourier component of EOM 



{kuj) 2 x k 



d(V) 



dx_ k 

Quantization of the action 
I=(p^)=n+^ 



Energy is phase average of H 
E(I) = (H) 
Classical frequency 

uj = dE{I)/dI 
Dynamical variable 
A(x(I,6),p(I,9)) 
Fourier coefficients 
A k (I)=f%A(I,e)exp(-ikO) 
Hamiltonian independent of angles 
Efc^o H k H-k = 
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V. FURTHER STUDY OF THE COMMUTATION RELATIONS 



A. One-dimension and WKB quantization rule 

In the following lines we present a derivation of the WKB quantization rule for bound 
states that is simpler than the one that we have previously published [2~0|. 

A diagonal element in a state n of the canonical commutation relation may be written 
(using reality conditions) 

oo 

2Xnn>iPnn> = 1- (5.1) 

n'=0 

Change n — > n" and sum this equation from to n. Divide the resulting double sum into two 
terms. The first, defined by taking both sums from to n, vanishes by antisymmetry (the 
matrix elements of x are symmetric under exchange of indices, those of p antisymmetric). 
We are left with the sum 

n oo 

E 2x n » n >ip n » n > =n + l. (5.2) 

n"=0 n'=n+l 



With our standard definition, 



x n ,- n (^(n + n')) = x nn >, (5.3) 



and a change of indices 

n" = n — u, n = n + k — u, (5.4) 

Q5.2D becomes 

OO fc— 1 J 

2 EE^( n +9 fc ~ v)ip k {n + -k - v) = n + 1. (5.5) 

Equation ( |5.5| ) is still exact and can used in quantum calculations. Expanding about the 
semiclassical value Xk(n) and keeping terms that contribute at most to order n and to order 
unity, we find 



E 

k=l 



2kx k {n)ip k (n) + (kx k {n)ip k (n)) + ... 
an 



n + l. (5.6) 



(To obtain the form of the second term on the left hand side of this equation requires a 
careful examination and grouping of terms that contribute.) The way to read this equation 
is to understand that the first term contains contributions of order n and smaller, whereas 
the second term is at most of order unity, and at the same time is half the derivative with 
respect to n of the first term. It follows that this second term has the value |. We thus 
derive the one-dimensional Cartesian WKB quantization rule 
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OO J 

2 ^ kx k (n)ip k (n) = n + -. (5.7) 
fc=i ^ 

A similar tour de force does not work in two or more dimensions. We omit the uninspiring 
details. The moral of the story is that we are "lucky" to have the EBK quantization rules 
for semiclassical quantization. Our aim when we started the investigation of this section was 
to discover an alternative semiclassical scheme based directly on the commutation relations. 
The conclusion, which we believe to be firm, is that this is not possible except in the one- 
dimensional case. 



B. Remarks on quantum aspects 

Consider a two-dimensional system and suppose it to possess only bound states. (In any 
event we use a notation with discrete labels only.) We can write a Heisenberg scheme that 



should be valid whether we deal with the "regular" spectrum or the "irregular" one |17 
We simply order the energy levels of a two-dimensional system with coordinates x, y, in a 
linear sequence, N = 0,1, ... . A calculation can be based on the equations 



{Ejfi — En) xnn' — —(F x , 



NN' 



(E N > — E N ) 2 y NN i = —(F y ) NN >, (5.8) 

yxffiv - e n )\x NN i\ 2 = i, 

N> 

J2(E N ,-E N )\y NN ,\ 2 = l, (5.9) 



A" 



(5.10) 



where En are the exact eigenvalues. These equations can be viewed in two ways, either 
as a set of sum rules to be satisfied by the exact solutions of the quantum mechanical 
problem, found by some other means such as matrix diagonalization in a basis, or else as the 
foundation for a computational scheme involving the solution of non-linear equations. For a 
well-defined quantum system (Hamiltonian bounded from below), we shall certainly be able 
to do a diagonalization and subsequently check the above equations. Since this should be 
possible both for the regular and for the irregular spectrum, this raises the hope that some 
version of the quantum matrix method can also be applied to the chaotic regime. 

toward this end, we note that by summing the commutation relations from to N, we can 
replace ( |5.9|) by the positive sums (this involves the same argument as in the one-dimensional 
case) 

oo N 

Y, E (En'-E n »)\x n »n'\ 2 = N + 1, 

N'=N+1 N"=0 
oo N 

Y E (E N ,-E N ,,)\y N ,, N ,\ 2 = N+l. (5.11) 

N'=N+l N"=0 

These relations are interesting because they guarantee the convergence of certain sums, and 
thus imply that the matrix elements cannot spread out too far as a function of energy 
differences. 
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VI. APPLICATION OF MATRIX MECHANICS AND COMPARISON WITH 

SEMICLASSICAL APPROXIMATION 



In this section we present illustrative applications of the matrix mechanics method for 
several simple models and compare the results with those of the semiclassical approximation 
as well as with the results of exact diagonalization. With a given semiclassical approxi- 
mation, as defined below, we shall associate the matrix method, that will be described as 
a sequence of approximations that should approach more and more closely to the exact 
result. Therefore we may anticipate that at a sufficiently high order of approximation its 
accuracy will exceed that of the semiclassical result. It is not, however, guaranteed to do 
this automatically for the simple reason that the quantum method that we shall describe 
does not have a true semiclassical limit, but is specifically designed to be most accurate for 
small quantum numbers. For this reason, we shall also describe an alternative method that 
starts with the EBK theory and specifically calculates quantum corrections to it. Numerical 
results for all three methods are presented at the end of this section. 

Imagine that we have a two-dimensional system described by a Hamiltonian 

H = ^ + ^ + ^(x 1 ) 2 + ^A 2 (a; 2 ) 2 + V fa * annaoSB (x\x 2 ), 

2 I \ |2 

= y + 2 + y anharmonic ( x ). (6.1) 

However, the methods described here are more general, and can be used with a system 
with many degrees of freedom. We shall describe these methods, in general, by referring to 
an n-dimensional system and give examples and show results for the specific cases n — 1 
and n = 2. For improved clarity, we have shifted the coordinate index from subscript to 
superscript. 

For the solution of the algebraic formulation of the semiclassical approximation we utilize 
an iterative method described by Percival and Pomphrey []I5,IS|. For the corresponding fully 



quantum calculations, we shall describe a quantum extension of this method. We have also 
carried out calculations using the Newton-Raphson method, but these will not be discussed 
here. 



A. Semi-Classical Iterative Method 

Following Percival and Pomphrey, this method utilizes the semi-classical equations of 
motion, given for Fourier components k (in terms of the above Hamiltonian): 

A-(k-^) 2 ]x k = F k (x), (6.2) 

where 

F(x) = _ W anharmonic ( x ), (6.3) 

and oj = (o4x,042, ■ ■ ■ ,oj n ) is the vector of frequencies of the system. When there is no an- 
harmonic driving force, uf = Aj. The vector k is an n-dimensional sequence of integers, 
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positive and negative: k= (ki, k 2 , ■ ■ ■ , k n ). It is also useful, here, to define the vectors 
1 = (1, 0, 0, . . . , 0), 2 = (0, 1, 0, . . . , 0), i = (0, . . . , 1, . . . , 0), where the "1" is in the i'th place 
and all other indices are zero. Thus, k • i = ki. In order to solve the equations that describe 
the semiclassical limit in algebraic terms, we must truncate the Fourier series by specifying 
a maximum integer vector |K|. The meaning of the absolute value is that we define this 
limit symmetrically with respect to the origin in k space. We set x\^ = for |k| > |K|. 

The procedure begins from the harmonic limit, where the only nonzero Fourier compo- 
nents are x l ±i . As an example, for a two-dimensional system, where we write x l = x and 
x 2 = y, the only nonzero components, initially, are x±i ) q and ?/o,±i- All others vanish. We 
shall refer to these as the driving components. The values of these initial nonzero compo- 



nents are computed from the action conditions, Eq. (|4.8|) , namely, x l ±i = Jli/2uii, where U 
is the action of the i'th degree of freedom, as given by Eq. ( |4.7| ). 

With these initial values, we can then calculate initial values for the Fourier components 
of the non-harmonic part of the force. Improvements to the frequencies u>i are then obtained 
from the equations of motion for the driving components, 

= A, - 3. (6-4) 
x\ 

With these new and improved frequencies, other Fourier components can be obtained from 
the remaining equations of motion. In particular, we compute 

x \ = kl > , for k ^ i. (6.5) 

k V-(k-cu) 2 T ^ 1 

This completes one cycle of iteration. We now return to the EBK quantum conditions, 
and using the calculated values of the frequencies and of the non-driving Fourier components 
we compute new and improved values for the driving components. Explicitly 



x\ = 

\ 



— [I i -5^A*(k-w)x k x_ k ]. (6.6) 



For the special case n 


= 2, 


we have 








£±1,0 = \ 


4- 


Z)(fci,jfc 2 ) ? £(±l,0) 


\ X k!,k 2 




ik\U)\ + k 2 UJ2) 




Vo,±i = \ 


ly- 


£(fci,fca)^(0,±l) ^ 2 




+ VkiM' 


[kiuj 2 + k 2 uj 2 ) 


2uj 2 



(6.7) 



(6.8) 

Once more we calculate the improved values of the frequencies and then of the remaining 
components and continue to loop until the energy converges. 



B. Quantum Iterative Method 

In this section, we discuss the modifications to the semiclassical perturbative method 
which are necessary in order to have a completely quantum calculation. The major compli- 
cation of the quantum approach is that in our version of the Heisenberg method, we must 
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consider the coupling of eigenstates both in the equations of motion and in the commutation 
relations. To deal with this problem in a practical way will necessitate the introduction of 
two cut-off parameters, compared to the one used in the semiclassical case. Consider a typ- 
ical matrix element, (m|0|n), of an operator O (in general a coordinate). If we have chosen 
an ordering of the integer vectors, we can suppose that n > m is defined and k = n — m is 
the analogue of a Fourier component. Since the first state always considered in the quantum 
method to be described is the ground state m = 0, K coincides with the maximum value 
of n that we include in the calculation. However, this specification does not completely 
define the calculation. The crudest way to define a closed scheme for fixed K is to allow 
only matrix elements with m = 0. But the equations of motion for the matrix elements 
(0|O|n) will necessarily bring in matrix elements with m 7^ 0. In order to have a closed 
scheme involving only the quantities specified, we must introduce a closure approximation 
that relates the unwanted matrix elements to those that belong to the allowed set. The 
crudest such approximation is to relate matrix elements by an equal displacement of both 
sets of quantum numbers, i. e., 

(m|0|n) = (O|0|n- m). (6.9) 

(A better approximation at this point is to use harmonic oscillator results for the large 
matrix elements that are non- vanishing in the uncoupled limit.) This crudest of closure 
approximations is called the zero band width (B = 0) approximation. 

To improve the scheme, we increase the band-width B. Thus in a B = 1 approximation, 
we treat without closure those values of m that comprise both the ground state and those 
states that are coupled to the ground state by harmonic oscillator matrix elements. By 
suitable extension we can define a scheme for arbitrary B, with the understanding that the 
maximum value of m, called M, cannot exceed K. A closure approximation, generalizing 
(p.9|), needs to be applied now only to relate matrix elements outside the expanded scheme 
to those in which the left hand state is M. As B increases, we find that the values obtained 
for low-lying states become more and more accurate, as well as insensitive to the closure 
approximation. The values of B and of K determine the number of equations of motion and 
commutation relations that have to be used to evaluate the set of included matrix elements 
and energy differences. This number is a small multiple of B times the number of equations 
for the semiclassical approximation associated with the same value of K. 

Turning to the detailed scheme, we begin from the harmonic limit. In this limit, all 
the energies are equally spaced. Thus, if we define quantum frequencies to be the energy 
differences 

{n\H\n} - (m\H\m) = E n - E, 
then the harmonic limit is defined by 

Here, the vector i is the same as was introduced above, namely, it is an n-dimensional vector 
whose entries are all zero except for the z'th one, whose value is unity. 
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(6.10) 



From the well-known solution for the harmonic oscillator, we also know the driving 
matrix elements x l n+i n . It is integral to our method that this result need not be put in 
from the outside but follows by application of the canonical commutation relations, given 
the result expressed by ( 16.111) . Actually the full statement is that for the harmonic case, as 
is well known, the simultaneous solution of the equations of motion and of the commutation 
relations yields the starting values of the quantum frequencies and of the driving matrix 
elements. In turn, we can calculate the lowest approximation to the force F* n+i . This 
completes a single cycle of the quantum iteration procedure. 

Next, the quantum frequencies are improved by using the generalization of Eq. (|(j. 1|): 



^n = A,-^N. (6.12) 

x n,n+i 



For further work, we write 



(6.13) 

a = ail + a 2 2 + ... + a n n. (6-14) 

It is straightforward to write this energy difference as a sum of fundamental energy differences 
or quantum frequencies in which just one of the quantum numbers changes by a unit. This 
is important in counting that the number of variables in the problem is determined by the 
equations of motion and a set of diagonal elements of the commutation relations. 

With the new and improved energy differences, it is possible to compute the perturbative 
quantum matrix elements from the equations of motion: 

<,n+a = "'^ fOT a ^ 1 ( 6 - 15 ) 

Lastly, we construct new "fundamental" matrix elements from the diagonal matrix ele- 
ments of the canonical commutation relations. Toward this end, we solve the n equations 
{% = l,...,n) 



i 

j..»"" ,,, - J 2 



W m j = - (6.16) 



for x l m+im . These equations replace the quantum conditions for the action variables, used 
in the semi-classical approximation, and, as already stated, play the same role as the latter 
did, in providing additional equations needed to have a determined set. 

We have thus identified the three distinct elements that are used in an iteration cycle: 
the commutation relations, the equations of motion for the "driving" matrix elements, and 
the equations of motion for the "small" matrix elements. 

C. Quantum corrections to the semiclassical approximation 

We start by reminding the reader that the semiclassical approximation to the energy, 
E{n) = (n\H\n) (in a one-dimensional notation), is obtained by expanding energy differences 
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and quantum matrix elements about certain values that are identified as semiclassical ones 
and dropping higher order terms. With the definitions 

u;(n)=E(n + ±)-E(n-±), (6.17) 



k = n"-n\ (6.18) 
1 

5< 



n=h n ' + n"), (6.19) 



we have to second order 



dxk{n) 1 2 d 2 Xk{n) 



(n'\x\n") = x k (n) + (n — n) — ; h -(n - n) — , (6.20) 

an 2 d z n 

E(n + k)- E(n) = ku{n) + + - (6.21) 

where the derivation of the second of these equations requires attention to the definition 

To obtain quantum corrections to the energy, we must retain and compute the values of 
the first and second derivative terms in (|6.20| ) and ( 6.21| ). These equations are expansions 
in (1/n) and we need terms up to second order because, as we have shown, the first order 
terms do not contribute to the energy. 

The calculation is straightforward. Reverting to an n-dimensional notation, we recall 
that the semiclassical quantization program is based on the equations 

(k-u) 2 xl = Fl, (6.22) 
fci(k • u)x{x 3 _ k = (n; + -Oi), (6.23) 

where F % k is a Fourier component of the ith component of the force and aij is a Maslov index. 
Let us write these equations succinctly as 

= 0, (6.24) 

where a is a label that enumerates in turn the equation set (|6.22| ) and ( |6.23| ), putting 
them in one to one correspondence with the Fourier components and with the frequencies. 
Correspondingly, z = {z a } = {x\,u>i}. 

We now treat these equations as if they are true for a continuous range of n and compute 
the first and second derivatives, leading to the equations, 

S a i, (6.25) 



r\ r\ '-'(I'll 

OZb OTli 

d^g d 2 Z b d 2 V a dz b dz c 

dzb driidrij dzhdz r drii dn 4 



(6.26) 



Here the column matrix s ai , consisting of all zeros and a single value of unity has its origin 
in the right hand side of Eq. ( |6.23| ). Both of the above sets of equations are linear and 
inhomogeneous and will have a solution provided the matrix 
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' 6 " 0* 



(6.27) 



is non-singular. (See the next section for further discussion of the uses of this matrix.) The 
solutions of these equations provides us with the necessary input for calculating the leading 
quantum corrections to the energy. 



D. Numerical illustrations 

We consider first the one- dimensional inverted quartic potential. 

H= -(p 2 + x 2 ) + bx\ (6.28) 

with b < 0. For this model as well as for the two dimensional model considered below, 
it is understood that these do not exist as quantum mechanical systems. Nevertheless 
diagonalization in a large - but not too large - basis will yield eigenvalues that we can treat 
as those of a bound system. Correspondingly the approximation methods described above 
will also yield values that we can take seriously and compare with the "exact" results. In 
the following and for the two dimensional model, we shall present a few examples of the 
calculations that we have done. 

In Table II we display the energies of the ground state and of two excited states calcu- 
lated for several values of b, both for the inverted quartic and the stable quartic oscillator. 
Four different calculations have been carried out at each point: E^g results from matrix 
diagonalization ("exact" value); E sc is the semiclassical result; E SC + AE is the semiclassical 
result with leading quantum corrections; E q is the value obtained from Heisenberg matrix 
mechanics. For the latter, the band width was chosen, when feasible, to give results in 
essential agreement with -Ediag- This required values of B = 3 — 5. Our general expecta- 
tion is that the semiclassical calculation with quantum corrections will be an improvement 
over the purely semiclassical calculation. This expectation is born out with the exception 
of the value for b = —.005, n = 10. Since the maximum of the potential for b < is at 
= 12.5 for b = —0.005, the n — 10 level gets sufficiently close to this maximum that 
the assumptions upon which our derivation of the semiclassical approximation and of the 
quantum corrections to it are based become suspect. It was to illustrate this point that 
results for the stable quartic potential were included. Here there is no sign of the difficulty 
encountered with the inverted potential. 

TABLE II. Representative values of the energy for states n of a one dimensional quartic os- 
cillator, for several values of the anharmonicity parameter b. -E^diag denotes the energy calculated 
by direct diagonalization, E sc is the semi-classical approximation, E sc + AE includes the lead- 
ing quantum correction to the semiclassical result, and E q is the quantum result obtained by the 
Heisenberg method. 

b n -Ediag E sc E sc + AE E q 

0.498489 0.499248 0.498500 0.498489 

-0.002 5 5.405378 5.406235 5.405521 5.405377 
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10 10.145887 10.146877 10.146244 10.145887 








0.496182 


0.498112 


0.496249 


0.496182 


-0.005 


5 


5.249313 


5.252083 


5.250752 


5.249314 




10 


9.484437 


9.489854 


9.493750 


9.484438 




n 
U 


u.oui4yu 


U.OUU/40 


n r:m r;nn 


n £m a nn 

u.oui4yu 


+0.002 


5 


5.588750 


5.588081 


5.588840 


5.588750 




10 


10.813669 


10.813058 


10.813814 


10.813668 







0.503687 


0.501862 


0.503746 


0.503687 


+0.005 


5 


5.712896 


5.711448 


5.713323 


5.712896 




10 


11.231303 


11.230091 


11.231885 


11.231304 



In Fig. 1 we compare, as a function of the coupling constant, semiclassical matrix elements 
with the two neighboring quantum matrix elements of which it should be the approximate 
average. (The notation makes it clear that these results are for the reference state n = 10.) 
The expectation is well satisfied for the largest matrix elements. We also see that except 
for the fundamental Fourier component, all other matrix elements are rising exponentially 
with increasing coupling constant, promising an eventual breakdown of the theory. 



26 



FIGURES 

V= 1/2 x 2 + bx 4 
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FIG. 1. Comparison of semiclassical matrix elements (Fourier components) for the inverted 
quartic oscillator and the reference state n = 10 with the quantum matrix elements that they most 
closely approximate, as a function of the anharmonicity parameter b. The square of the value is 
plotted with the semiclassical result shown as as a dashed line and the associated quantum values 
as full lines. 

We turn next to the two dimensional generalized Henon-Heiles model, 

H = l -{pl + Ax 2 ) + \{ V \ + /V) + V(x, y), (6.29) 
V(x,y) = by(x 2 + b'y 2 ), (6.30) 

where we choose A = 1.69 and \i = .49, in order to stay away from low-order resonances. We 
shall consider a range of values of b and b' in order to expose different physical situations. 

In Table III we display a limited number of results analogous to those presented in Table 
II. Here the failure of the quantum corrections to the semiclassical result is flagrant for the 
first choice of interaction parameters, which can be ascribed to the smallness of the height 
of the lowest barrier. This barrier is raised by the second choice of parameters, and this 
seems to cure the difficulty. Again, the accuracy of the non-linear calculations carried out 
for advertised form of Heisenberg matrix mechanics is to be noted. 
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TABLE III. Energies of states labeled by n\,ri2 for the Henon-Heiles model denned in the 
text, as a function of the anharmonicity parameter b, for a fixed value of b' . E^ ag denotes the 
energy calculated by direct diagonalization, E sc is the semi-classical approximation, E sc + AE the 
quantum-corrected semiclassical approximation, and E q is the quantum result obtained by the 
Heisenberg method. 



b b' 


ni n 2 


Ediag 


E sc 


E sc + AE 


E g 







0.999628 


0.999627 


0.999757 


0.999628 


-0.04 0.01 


2 2 


4.989952 


4.990620 


4.991267 


4.989952 




4 4 


8.967166 


8.969413 


8.969985 


8.967164 







0.998504 


0.998501 


0.998817 


0.998504 


-0.08 0.01 


2 2 


4.958779 


4.961547 


4.961127 


4.958780 




4 4 


8.862238 


8.871849 


8.862532 


8.862098 
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In Fig. 2 and Fig. 3 we compare semiclassical and quantum matrix elements for the 
coordinates x and y, respectively, in analogy with what is shown in the one dimensional case. 
The same general remarks apply here as for that case. Note that the range of parameters 
does not include those used in Table III. 
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FIG. 2. For the potential shown at the top of the figure, comparison for the state {n\ = 712 = 4) 
of quantum values and of the associated semiclassical approximation for selected matrix elements of 
the coordinate x. The full lines represent the quantum results and the dashed ones the semiclassical 
values. 
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FIG. 3. The same comparison as in Fig. 2, for the coordinate y. 

VII. ADDITIONAL OBSERVATIONS AND SUGGESTIONS CONCERNING THE 
SEMICLASSICAL AND CLASSICAL ANALYSES 



In this section, we raise a number of issues related to the semiclassical analysis that it 
may to profitable to explore. Since some of the suggestions put forward in this section are 
speculative, they may also turn out to be worthless. 



A. Fixed action quantization program 

We first consider other possible applications of linear variation about the semiclassical 
equations, as embodied in Eq. ( |6.25| ) . These equations are a special case of the equations 

S a ,b$Z b = S a , (7.1) 

which has a number of possible applications, depending on the choice of the source vector 
s. We discuss briefly several such applications. 

Al) There is first of all the basic algorithm of the Newton-Raphson method. In this 
instance, we choose s a = —^y, corresponding to a vth approximation to z a , and the 
matrix S is also known in this approximation. The solution to ([7. 1|) determines the next 
approximation to z a , namely 
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(7.2) 



A2) If the above solution converges, we thereby define a "stability matrix", S. As we 
have already seen, it is this quantity that occurs when we differentiate Eqs. fl6.22| ) and ( |6.23| ) 
with respect to rij, leading to Eqs. (|6.25|) and ( |6.26|) . In addition to the application already 
described in the previous section, we can think of at several other applications of these 
equations. We consider only ( |6.25|) , but ( |6.26|) could be brought in. For example, they allow 
us to compute the energy of neighboring states, according to the formula (two-dimensional 
example), 



E(ni + d\, n 2 + d 2 ) = E(n 1 ,n 2 ) + u){d{ + ^tUijdidj, 



(7.3) 



where = (didi/drij). Work found in the literature, as far as we are aware, uses only the 
linear approximation to this result. As a second application, when we want to calculate 
solutions for a fixed Hamiltonian and for values of the action neighboring to those for which 
solutions are already known, the solutions of the linearized equations can be used to obtain 
improved starting values for the Newton-Raphson iteration. 



A3) To obtain solutions for a neighboring Hamiltonian, one can again use a form of ( pM| ), 
with yet another driving term easily derived from the structure of the original nonlinear 
equations. 

An important question that should be susceptible to study by the linearized formalism 
is the relation of an instability of a solution of the nonlinear equations to the eigenvalues of 
the matrix S. 



B. Solutions at fixed frequency; applications 



We now restrict our attention to the solutions of the equations of motion ( |6.22j ) for fixed 
frequency, i. e., we study the purely classical problem, setting aside for the moment the 
question of how to adjoin a quantization feature. Let us imagine that we are interested 
in obtaining the Fourier coefficients as functions of the frequency values. There arises the 
practical problem of how to choose a sensible grid of frequency values. This problem can be 
solved presumably by perturbation theory for small coupling, and we can use a form of the 
linearized analysis described above to show us how to change frequencies locally. As a first 
approximation to quantization - if only energies are of interest - we can compute values of 
the action from the Fourier coefficients and the frequencies and write 

Ji = n + ^ (7.4) 

where we have used the tilde to indicate a non-quantized value, so that rj is not an integer. 
We can then extrapolate to the nearest quantized values of the energy by using Eq. ( [7.3| ). If 
we want also the Fourier amplitudes for the quantized invariant tori, we have to successively 
modify the frequencies until they assume values associated with quantized actions, for this 
purpose we might use the formula, 
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duo- d H 

^ = dt 6J 3 = dJ^^- r i ] - (7 - 5) 

The most practical way of calculating the required second derivatives would probably be 
from a grid of energy values. 

We finally note another possibly interesting way of utilizing a grid of fixed frequency 
solutions to quantize a system. Consider, as a simple example, a two dimensional non- 
resonant system. Write the classical Hamiltonian in the normal form 

H = Y,a mim2 JrJ? 2 . (7.6) 

It would appear at first sight that the number of points on our grid of solution values will 
determine the number of terms that we can use in this equation, which yields a set of lin- 
ear inhomogeneous values for the coefficients a mim2 . We would then quantize the resulting 
Hamiltonian by using the EBK quantization conditions or some variant, following the dis- 
cussion of Ref. If we consider the usual procedure for constructing a form such as (|7.6|), 



however, we would guess that our previous remarks are much too naive. One's ability to 
carry the expansion to a higher and higher order is contingent upon obtaining a perturba- 
tion expansion to the appropriate order. To make contact with methods based on Fourier 
series, we must note that an analysis can be carried out which informs us which Fourier 
components must be included to guarantee equivalence to a perturbation expansion up to 
a prescribed order. The solution of a nonlinear scheme including only these components, 
at the same time that it contains a selective (and uncontrolled) summation of higher order 
terms, is at least perturbatively correct to some controlled order. It is this latter order which 
would determine how many terms are allowed in the expansion ( |7.6| ) • To our knowledge this 
relationship between the methods based on Fourier expansion and those based on normal 
forms has not been considered previously. 



C. Direct use of the variational principle 

Can one use the variational principle directly to simply solution of the nonlinear equa- 
tions? Let (L)^ be the vth approximation to the average Lagrangian. Expanding to first 
order, 

, I) = ,I) M + f%,A (7.7) 
dx it k OuJi 

Consider first the fixed frequency case and choose 

5x i)k = c-H— • (7-8) 



5xi 



i . k 



This is useful provided the average Lagrangian is truly an extremum. It would then appear 
that a suitable choice of the constant c u both as to magnitude and sign would move the 
average Lagrangian to its extreme value. The reason that we must allow the constant to 
depend on the order of approximation is that we want to guarantee that the retained first 
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order correction is larger than the omitted second-order terms. This means that we must 
start out with conservatively small values of the constant and let it increase toward unity 
as (if) we approach convergence. 

If we carry out a calculation at fixed action, we have to consider the equation, 



8{L) 



4 U) = Ji- (7.9) 



8uji 

Comparison with ( |7.7| ) suggests that we choose 

5ujj = bjJj (no sum). (7-10) 



From this assumption, the assumption of fixed J, and ( |7.8| ) we can derive a pair of linear 
equations for bj, namely 

^ dJi dJi 8{L) 



j 



dujj dxj t k 8xj t k 



Here indices indicating the order of approximation are omitted. In the present case, there 
is no guarantee that the two first order correction terms in (|7.7|) are of the same sign, and 
thus the way that convergence may be achieved, if at all, is somewhat more problematical. 



VIII. CONCLUDING REMARKS 

In this paper we have presented a new view of the transition from Heisenberg matrix me- 
chanics to the theory of invariant tori. We have suggested a number of possible applications 
of the ideas that were presented mainly in Sees. II, IV, and VII. We have applied the ideas 
developed in Sees. II and IV to several standard models. The most important feature that 
has emerged from these applications is that the Heisenberg mechanics can be developed into 
a quantum calculus that is only modestly more complicated than the semiclassical calcula- 
tions associated with the theory of invariant tori. This can be done in two ways, either by 
starting with the semiclassical approximation and building a correction scheme about it, or 
else by constructing a fully quantum scheme starting from the ground state. 

If one were now to ask for the most important next step that one could take with the 
Heisenberg methods, an excellent candidate for an answer would be to produce solutions 
for a globally chaotic system such as the one studied by Martens et al [|||. It would also 



be worthwhile to revisit some old ground. As an example, we might restudy the edges 
of the regions, as a function of coupling strength and quantized actions, beyond which 
converged solutions of our equations cannot be found. Though some of this was done in 
the present work, the issue is somewhat muddied for the model chosen, since it does not, 
strictly speaking, possess a Hilbert space. There is evidence, based on calculations for the 
standard mapping [P7| , P%fl , that this failure represents an independent method for studying 
the disappearance of invariant tori. 

Another idea that might be reexamined is that of approximate tori. Reinhardt and as- 
sociates have championed this idea, that permits EBK quantization to be applied past the 
point of breakdown of the associated invariant torus. This idea is suggested both by their 
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work on the quantization of normal-form Hamiltonians [^9|,[3(J and on adiabatic switching 
23fl . It reappears in the work of Martens and Ezra [32,331 in that trajectories that appear 



to be associated with locally chaotic regions show enough well-defined Fourier components 
so that quantized actions can be calculated. Independently, we realized that this phenom- 
ena would manifest itself in our work as follows: The non-linear equations underlying our 
approach yield a solution when only a small number of Fourier components are retained, 
but the solution blows up when an attempt to add components is made. By contrast, when 
invariant tori exist, this is signaled by insensitivity of the solution to the addition of Fourier 
components beyond a fixed number. The finite Fourier sums describe approximate invariant 
tori in the same general sense, though in a different approximation, as the approximate 
normal-form Hamiltonians. 

The transition from the quantum to the classical domain by the methods of this paper 
presents a new aspect of the study of the consequences of the order in which the two limits 
h — > 0, t — > oo are taken For the order studied in this paper, in which the time limit is 
taken first, it is quite impossible to strictly reach the regime of multiply periodic motion, as 
has already been pointed out in the body of our work. It may be of interest to to undertake 
a further study of the equations that can be obtained in this limit. There may also be some 
connection of these ideas with the idea of approximate tori. 
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APPENDIX A: PROOF THAT SOLVING THE EQUATIONS OF MOTION 
DIAGONALIZES THE HAMILTONIAN 

We prove directly that a solution of the equations of motion and of the commutation 
relations guarantees that the Hamiltonian has been rendered diagonal. From the equations 
of motion, we can derive, in a way which is obvious from the right hand side of what follows, 
the equation 



J2{(Ei-E m ) 2 + {E, 



E x ) 2 }x ml x lm = (xV + V'x) 



mn 



2{xV) 



(Al) 




mn 



yXffi — E m )(E n — Ei) XmlXl 



(A2) 



Equations ( |A1| ) and ( |2.5| ) can be combined in several useful ways: Thus by subtracting twice 
(p|) from Q) we find 

E n - E m ) 2 {x 2 ) mn = -2{p 2 ) mn + 2{V'x) mn , (A3) 



which is the matrix element of the equation 



^[[x 2 ,H],H] = -p 2 + V'x). 



(A4) 



On the other hand by adding twice ( p.5|) and (|A1|) , we find the results 



£)(2£, -E m - E n ) 2 x ml x ln = 2{p 2 ) mn + 2{V'x) 



(A5) 



Equation (|A~5|) will be used to prove that H is diagonal in conjunction with another relation 
which is a further consequence of the equations of motion, namely, 
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{E n - E m ) £(2£, -E m - E n fx ml x ln = -2i{V'p + P V') mn + 2{E n - E m ){V'x) mn . (A6) 
i 

Since 

i(V'p + P V') mn = 2[V, H] mn = 2{E n - E m )V mn , (A7) 

we have finally, for m ^ n, in place of ( |A5|) , 

^2(2E t -E m - E n fx ml x lm = -AV mn + 2{V'x) mn . (A8) 
l 

Comparing ( |A4| ) with (|A8[), we conclude that 

2{p )rrm 2(V x) m n 4V^j n -|- 2(V x) mrl , (A9) 

or 

Hmn = 7,(,P 2 )mn + V mn = 0, 771 ^ 71. (A10) 
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